function kap_return = KAPPA(a,b,c,x0)
kap_return=zeros(size(x0));
 y_prime = @(x) 3*a*x.^2 + 2*b*x + c;  
 y_double_prime = @(x) 6*a*x + 2*b;      
    for i = 1:length(x0)  
        x = x0(i);  
        y_prime_value = y_prime(x);  
        y_double_prime_value = y_double_prime(x);  
        kap_return(i) = abs(y_double_prime_value) / ((1 + y_prime_value^2)^(3/2));  
    end  
end  
